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Although it is well known that there are mutation hot spots in mtDNA, whether there are damage hot spots 
remain elusive. In this study, the regional DNA damage of mitochondrial genome after ionizing radiation 
was determined by real-time quantitative PCR. The mtDNA damage level was found to be dose-dependent 
and regional unequal. The control region was the most susceptible region to oxidative damage. GGG, as an 
typical hole trap during charge transport, was found to be disproportionally enriched in the control region. 
A total of 107 vertebrate mitochondrial genomes were then analyzed to testify whether the GGG enrichment 
in control region was evolutionary conserved. Surprisingly, the triple G enrichment can be observed in most 
of the homeo thermal animals, while the majority of hetero thermic animals showed no triple G enrichment. 
These results indicated that the triple G enrichment in control region was related to the mitochondrial 
metabolism during evolution. 

The human mitochondrial genome is a small (16569-bp) closed- circular, duplex molecule that is present at a 
multiple number per cell. MtDNA contains 37 genes including 13 structural genes for the components of the 
respiratory- chain enzyme complexes, 2 ribosomal RNAs, and a complete set of 22 tRNAs required for the 
translation of mtDNA-encoded mRNAs 1 . MtDNA is believed to be more susceptible to DNA damage and 
consequently acquires mutations at a higher rate than does nuclear DNA. The differences in its susceptibility 
could be caused by lack of protective histones, limited DNA repair capacity, and high level of reactive oxygen 
species produced during oxidative phosphorylation 2 " 4 . It is estimated that the mutational rate of mtDNA is at least 
10 times higher that of nuclear DNA. High polymorphism is characteristic of the mitochondrial genome. 
Transitions C-T and G-A occur at the rate of 2*10" 7 per site per year in mammals 5 . On the other hand, GC- 
rich mtDNA sequences and especially palindromic regions are highly sensitive to mutagenesis 6 ' 7 . Especially, there 
are positions and short sequences with an extremely high mutation rate 8,9 . One of the direct causes of mutations is 
the specific DNA damage at certain positions. As the mutation "hot spots" are widely indentified in mitochon- 
drial genome, it is of interest to figure out whether there are mtDNA damage "hot spots" correspondingly. 

MtDNA damage can be measured by Southern Blotting, 8-oxoG damage, or a comprehensive scanning of the 
mitochondrial genome by RFLP or TTGE analyses. However, these approaches are labor intensive and require 
large amounts of DNA. To overcome the difficulties in mtDNA damage evaluation, Santos et.al. established a 
long-run quantitative PCR method that could permitted the amplification of the whole mtDNA genome 10 . The 
rationale is that any lesion on DNA strands will interfere the PCR amplification. On that basis, real-time QPCR 
was developed for better quantification of mtDNA damage 11 . The semi-long run real-time PCR approach 
allowing analysis of differential DNA damage in the mitochondrial genome 12 . Although region -specific 
mtDNA damage has been investigated by Rothfuss et.al. 12 , only four regions in mtDNA genome have been 
analyzed and the whole picture of regional mtDNA damage is still lacking. In addition to the reports of 
Rothfuss et.al., Merino et.al. found that damage and mutations are preferentially funneled to the conserved 
sequence block II within the control region of the mitochondrial genome through DNA charge transport 13 " 15 . 
Charge migration in DNA has been a long lasting interest of the radiation community due to its relevance for the 
mechanisms of DNA oxidative damage 16 . The relative simplicity of mtDNA due to its prokaryotic origin provided 
favorable conditions to investigate the mechanism involved in the regional mtDNA damage. 

Here, we developed a semi-long run real-time PCR method for the accurate quantification of mtDNA damage 
induced by ionizing radiation (IR). For the first time the regional mtDNA damage covered the whole mitochon- 
drial genome was investigated. The mtDNA damage was found to be dose-dependent and regional unequal. The 
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control region was found to be the most susceptible region to oxid- 
ative damage. A specific nucleotide motif, GGG was disproportion- 
ally enriched in the control region. A total of 107 mitochondria 
genome were then analyzed to testify whether the GGG enrichment 
in control region was genetic conserved. 

Taken together, the first mitochondrial genomic mapping of DNA 
damage induced by IR was reported. Experimental data and bioin- 
formatics analysis provided data supporting that charge transport 
play an important role in the radiation-induced region -specific mito- 
chondria DNA damage. The triple G enrichment in control region of 
homeothermal animals provided evidence of mtDNA evolution in 
vertebrates. 

Results 

The D-loop region is most vulnerable to radiation induced 
mtDNA damage. To investigate the regional mtDNA damage 
induced by ionizing radiation, 11 pairs of primers were used to 
amplify the whole mitochondrial genome. Each of the amplicons 
was about 1.8 kb, with one exception that was 1077 bp (region 
13902-14979). The results of mtDNA damage by SLR rt-PCR 
analysis were shown in Figure 1. Dose-dependent increase of 
mtDNA damage can be observed in all the regions analyzed. The 
SLR rt-PCR analysis revealed a similar lesion rate in all the encoding 
regions from position 1196 to 14979. Initial mtDNA damage in 
MCF-7 cells treated with 2 Gy X-ray irradiation ranged from 0.87 
lesions per lOkb DNA in the region 4487-6286 to 1.51 lesions per 
10 kb DNA in the region 7266-9077. A steady increase of mtDNA 
lesions was observed with the increase dosage of X-ray irradiation 
from 2-8 Gy. Strikingly, there were approximately 3.8-fold and 2.5- 
fold increase in region 14898-151 and 16488-1677, respectively with 
2 Gy irradiation. The significant elevation of mtDNA lesions in these 
two regions can also be observed in MCF-7 cells with 4 Gy and 8 Gy 
irradiation, although the fold increase was alleviated. Interestingly, 



both region 14898-151 and 16488-1677 overlapped with non-coding 
region of mtDNA (16024-576), while all of the less-damaged areas 
(1196-14979) were within the coding region of mtDNA. 

To testify whether the regional mtDNA damage was cell line 
dependent, HeLa and HepG2 cells were employed. Similar results 
of initial mtDNA damage in each region was observed in HeLa and 
HepG2 cells after 4 Gy irradiation (data not shown), indicating that 
the regional mtDNA damage was not cell line specific. 

To testify whether there was regional mtDNA damage after irra- 
diation, as shown by SLR rt-PCR, a classic mtDNA damage assay, 
Long-run quantitative PCR was employed. As shown in Fig. 1. B, the 
mtDNA lesions of the 9.4 kb fragment harboring the D-loop region 
were much higher than the 8.4 kb amplicon that containing no D- 
loop region after 4 Gy irradiation. These results confirmed that there 
were preferentially damaged regions in mitochondrial genome after 
ionizing radiation. 

The triple G was enriched in D-loop region. The complete map of 
mtDNA damage after oxidative insult provided evidence supporting 
that D-loop played a role of sacrifice during oxidative stress. So why 
most of the oxidative burden was loaded on the D-loop? It has been 
long recognized that radiation-induced DNA damage are not 
randomized 1718 . Due to the different oxidation potential of DNA 
bases, guanine was thought to be the most easily oxidized base 
among all four bases of DNA 19 . It was therefore reasonable to 
speculate that region enriched with GC content could be more 
vulnerable to oxidative damage. However, both the severely 
damaged regions (14898-151 and 16488-1677) showed no overall 
biased GC enrichment (48% and 45.4%, respectively). Another 
explanation was that the non- randomness for DNA damage 
induction could be dependent on the primary structure of the 
adjacent DNA sequences 20 . If so, specific nucleotide motifs, rather 
than individual nucleotide positions, are targets of oxidative insult. 
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Figure 1 | Radiation-induced dose-dependent and regional specific mtDNA damage. (A) Quantification of mtDNA damage per 10 kb DNA by SLR rt- 
PCR amplification of total DNA isolated from MCF-7 cells exposed to 2-8 Gy X-ray irradiation showing an steadily increasing mtDNA damage with 
increased dosage in all tested mtDNA regions. Error bars designate standard deviation (at least three independent experiments). Region 14898-151 and 
16488-1677 exhibited significant increased mitochondrial damage, overlapping the control region of human mitochondrial genome. (B) Quantitative 
measurement of mitochondrial DNA damage using Long-run quantitative PCR. Representation of the DNA lesions of the 8.9 kb and 9.4 kb 
mitochondrial fragments. Lesion frequencies of treated samples were calculated per amplicon size and expressed per 10 kb of mitochondrial genome. 
Error bars represent the SD, each done in at least triplicate. ^Statistical significant at p < 0.05; ^^Statistical significant at p < 0.01. 
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To investigate whether there were specific nucleotide motifs selec- 
tively accumulated in the most damaged region of mtDNA, compu- 
tational scanning of the motif combinations in mtDNA was carried 
out. Initially, the complete 16569 bp mtDNA was divided by 
1841 bpX9. The scanning began at the position 14898 according 
to the SLR rt-PCR results. The outputs showed that while most of 
the three code combination of nucleotides were distributed ran- 
domly along the mtDNA circle (Supplement materials), the GGG 
motif was unevenly accumulated in region 14898-169 (Fig. 2A). The 
enrichment of G- stretch in this region can also be detected by the 4-6 
codons scanning (data not shown). Since D-loop has been proposed 
to be the most vulnerable region to oxidative damage in mtDNA, 
scanning began at position 16024, which coincided with the begin- 
ning of D-loop, was carried out. As shown in Fig. 2B, significant 
increased number of triple G in D-loop region was detected. 

Triple G enrichment in control region was restricted in homeo- 
thermal animals. If charge transport does play an important role in 
mitochondrial biogenesis, it would not be a special case for human. 
To further confirm the biological importance of charge transport in 
mtDNA, phylogenetic analyses combined with the screen of triple G 
enrichment in control region were carried out. 107 vertebrate species 
were analyzed and the result was illustrated in figure 3. Strikingly, 
most homeothermal animals (30/31) tested here showed triple G 
enrichment in control region, which were ranging from 96.41%- 
473.09%. The average abundance ratio for homeothermal animals 
was 21 1.95%. while most of the heterothermic animals (71/76) showed 
lower triple G distribution in control region (less than 100%), as mea- 
sured by triple G abundance ratio (Fig. 3). The average abundance 
ratio for heterothermic animals was 56.30%. 

Discussion 

The mitochondrial genome is known to be more fragile than nuclear 
DNA under ROS attack. Moreover, there has been speculation about 
a variable vulnerability to oxidative stress along the mitochondrial 
genome. Rothfuss et.al. first reported that D-loop was more prone to 
ROS-derived damage than other three region along mtDNA 12 . In line 
with their study, our current results showed that control is the one 
and only hot spot along mtDNA, as the complete mapping of re- 
gional mtDNA damage was provided. The control region is the most 
polymorphic region of mtDNA genome, it is therefore not surprise 
that initial damage inflicted by IR was most extensive in this area. 
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Since mtDNA control region play a key regulatory role of mtDNA 
transcription and replication, the preferential funneling of damage to 
control region could served as a passive checkpoint for mtDNA 
dynamic control. 

In radiation society, it has been a decades-long mystery that radi- 
ation-induced DNA fragmentation was not randomized. Nuclear 
DNA was too complex to investigate, but mtDNA may provide a 
more simpler model to investigate the mechanism involved. In our 
study we found that GGG was selectively enriched in control region. 
Triple G acted as the hole trap during charge transport (CT) 21,22 . 
Accordingly, oxidative damage could be funneled to the triple G site 
via charge transport. Our current findings underscored the import- 
ance of charge transport on mtDNA damage and mutations, as pro- 
posed by Merino et.al. 1415 ' 23 . The preferential funneling of damage to 
the control region via CT may provide a checkpoint for mtDNA 
replication and transcription, since D-loop regulated the replication 
and transcription of mtDNA. In addition, with the positive charge 
has a high probability of migrating towards the non-encoding area of 
mtDNA, the gene essentially for the electron transport chain was 
thus protected against the mutations. 

The role of selection in mtDNA evolution attracted little attention 
until recently 24-26 . Current results pointed to the importance of 
charge transport in mtDNA during evolution progress. Homeother- 
mal animals need to sustain high metabolic rate for the maintenance 
of body temperature (36-38°C for most mammals and 39-42°C for 
most birds). Mitochondrial oxidative phosphorylation (OXPHOS) 
generates not ATP but also heat, which is essentially important for 
the maintenance of body temperature in homeothermal animals. 
Shen et.al. suggested that strongly locomotive birds would have 
experienced stronger evolutionary constraints to eliminate deleteri- 
ous mutations in mtDNA to maintain efficient energy metabolism, as 
compared to weak locomotive birds 27 . If so, homeothermal animals, 
which require more energy metabolism to sustain body temperature 
than heterothermic animals, should have also subject to stronger 
evolutionary constraints to minimize the deleterious mutations. 
The funneling of mtDNA damage to the non- encoding region via 
charge transport, thus provide additional weapon to defend the oxid- 
ative stress during energy metabolism in homeothermal animals. It is 
interesting to note that while triple G enrichment in the control 
region of mitochondrial genome was discovered in current study, 
triple G enrichment was obviously enriched at the end of linear 
chromosomes, called telomeres. Telomere repeat, TTAGGG, which 
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Figure 2 | Computational scanning of the triple G across the mitochondrial genome. The 16569 bp mtDNA was divided into 9 sections, each contained 
an equal number of 1841 bp fragment. The scanning began at position 14898 according to the SLR rt-PCR data (A) or began at position 16024 of the 
beginning of control region (B). The scanning program was design by C language. 
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Figure 3 | Computational score of the abundance ratio of triple G in D- 
loop region of 107 vertebrates. The scanning program was design by C 
language. The phylogenetic tree of vertebrates were constructed based on 
107 complete vertebrate mitochondrial genome available in GenBank. 



was evolutionary conserved in vertebrates 28 , was also susceptible to 
oxidative stress 29 . Further investigation is needed for better under- 
standing of the effect of charge transport on genome evolution. 

Methods 

Cell culture and irradiation treatment. The human breast cancer cell line, MCF-7, 
liver hepatocellular cells, HepG2 and cervix carcinoma cells, HeLa were purchased 
from the American Type Culture Collection. Cells were maintained in Dulbecco's 
Modified Eagle's Medium (DMEM; Gibco) supplemented with 10% fetal bovine 
serum. Cells were cultured in 5% C0 2 in humidified air at 37°C. 

X-ray were generated with an X-ray machine (FAXITRON RX650, USA) operated 
at 130 keV. An exposure-rate meter (AE-1321 M, Applied Engineering Inc, Japan) 
was used for the dosimetry. The dose rates were 1.3 Gy/min. Cells in exponential 
growth were irradiated at room temperature, with non-irradiated culture cells 
(control) which were handled in parallel with the irradiated samples. 

DNA isolation and purification. Samples after irradiation were immediately 
subjected to DNA isolation. Total DNA was purified using DNA Blood and Tissue Kit 
(Tiangen, China) from cells with/without ionizing radiation, respectively, and DNA 
quantity and purity was determined by spectrometric analysis. The isolated DNA 
showed a high purity (A260/A280 > 1.8) and was stored at 4°C according to standard 
procedures. 

Semi-long run rt-PCR. The SLR rt-PCR amplifications were conducted in the FTC- 
3000 qPCR system (FUNGLYN, Canada), and the amplification was monitored and 
analyzed by measuring the intercalation of the fluorescent dye to double- stranded 
DNA supplied by the MightyAMP ver.2 kit (Takara, Japan). The mtDNA regions, 
primers and their PCR efficiency for rt-PCR are displayed in Table 1. 

The PCR conditions for the different fragments were optimized to achieve similar 
amplification efficiencies required to compare different amplicons. The product 
specificity was monitored by melting curve analysis and sequenced by Sangon, 
Shanghai for identification (data not shown). The reaction mix (total volume V = 20 
|il) consists of 1 X MightyAMP real time mix (Takara, Japan), 500 nM each forward 
and reverse primer and the equivalent quantities of template DNA (10 ng of total 
DNA). The cycling conditions include a pre-incubation phase of 10 min at 95°C 
followed by 40 cycles of 10 s 95°C, 10 s 55°C, and 60 s 72°C Each sample was assayed 
in quadruplicate, fluorescence was continuously monitored versus cycle numbers and 
threshold cycle (Ct) values were calculated automatically. 

Data analysis is based on the measurement of the Ct. Isolated total DNA from 
untreated sample was taken as reference. For each of the 1 1 mtDNA regions the 
difference in the Ct ACt (IR/non-IR) was used as a measure of the relative mtDNA 
damage with the 2" ACt method in correlation to the amplification size of the long 
fragment 19 . The Ct of IR-damaged mtDNA versus the Ct of non-irradiated mtDNA 
was directly employed to represent the level of mtDNA damage in each region. The 
short segment of each region was not employed for reference due to the fact that even 
the amplification of short amplicon could be affected by IR (Figure 4). The DNA 
damage was calculated as lesion per 10 kb DNA of each mtDNA region by including 
the size of the respective long fragment and displayed as average of at least three 
independent experiments. 

Lesion rate[Lesion per 10 kb DNA] = [l-2~ Along ) x 
I0000[bp] / size of long fragment[bp] 



Long-run quantitative PCR for mtDNA damage assay, long-run quantitative PCR 
method was conducted to validate the results generated by the SLR re- PCR method. 
Two amplicons were employed, the 8.9 kb amplicons were from Santos and 
colleagues 10 , and additional 9260 bp amplicon was used for confirmation of region 
specific mtDNA damage. The primer information was listed in table 1. The 
experiment procedures were carried out according to the protocol of Santos et.al. 10 . 

Computational analysis of motif distribution in human mtDNA. A computer 
program was designed for the scanning of sequence motifs distribution on mtDNA. 
The input parameter included the length of motifs and the fragment size to divide the 
whole mtDNA. The output will be counting number of each motif in each fragment. 

Computational scanning of triple G abundance ratio in D-loop region of 
vertebrates. The GGG abundance ratio was determined by the ratio of the number of 
triple G of D-loop to the total number of triple G in the mitochondrial genome. The 
scanning program was design by C language. A total of 107 vertebrate complete 
mtDNA sequences were analyzed. 
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Figure 4 | Reduced amplification of different amplicons after ionizing 
radiation. Analysis was performed immediately after exposure to 
radiation. IR/non-IR ratios were calculated by dividing the IR-damaged 
mtDNA signal for amplicon by the untreated control signal and expressing 
the ratio as a percentage of the untreated control set at 100%. Error bars 
represent the SD, each done in at least triplicate. * Statistical significant at 
p < 0.05. 



Phylogenetic analyses. A total of 107 vertebrate mtDNA sequences were randomly 
chosen from NCBI nucleotide database, which fitted for both triple- G abundance 
ratio scoring and phylogenetic analysis. Sequences were aligned first using the 
program Clustal_X 30 and refined manually. Phylogenetic analysis was done using the 
MEGA4.0 31 . Distance matrices were constructed using the Kimura two-parameter 
model and trees generated by the neighbor-joining (NJ) distance method. The 
bootstrap values were estimated from 1,000 replications. 

Statistical analysis. Statistical analysis was performed on the means of the data 
obtained from at least three independent experiments. Data are presented as 
means ±SD. Student's t-test program in Microsoft Excel was used to detect statistical 
significance. p<0.05 was considered to be statistically significant. 
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